Spectra of primordial fluctuations in two-perfect-fluid regular bounces 
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We introduce analytic solutions for a class of two components bouncing models, where the bounce 
is triggered by a negative energy density perfect fluid. The equation of state of the two components 
are constant in time, but otherwise unrelated. By numerically integrating regular equations for 
scalar cosmological perturbations, we find that the (would be) growing mode of the Newtonian 
potential before the bounce never matches with the the growing mode in the expanding stage. For 
the particular case of a negative energy density component with a stiff equation of state we give a 
detailed analytic study, which is in complete agreement with the numerical results. We also perform 
analytic and numerical calculations for long wavelength tensor perturbations, obtaining that, in 
most cases of interest, the tensor spectral index is independent of the negative energy fluid and 
given by the spectral index of the growing mode in the contracting stage. We compare our results 
with previous investigations in the literature. 

PACS numbers: PACS numbers: 98.80.-k, 98.80.Cq, 98.80.Es 
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INTRODUCTION 

Cosmological models with a bounce [li - a contraction 
which reverses into an expansion - may solve the hori- 
zon problem |2:1 in a non-inflationary way, i.e. as the 
Pre Big-Bang [3] and Ekpyrotic scenarios. In order to 
make bouncing models real competitors of inflation, or 
at least complementary to it, while addressing the singu- 
larity problem unavoidably present in such models, the 
spectrum of density perturbations which results from the 
bounce should be understood as it is in inflationary the- 
ories. Unfortunately, the physics of cosmological pertur- 
bations during a bounce is much more subtle, because of 
the reversal of growing and decaying modes before and 
after the bounce, so that even though the bounce dura- 
tion itself may be very short, usual matching conditions 
[^, @| should be used with particular care and verified. 

These subtleties have generated many works on the 
subject, in particular after the proposal of the Ekpyrotic 
scenario Q , which is based on a very slow contraction and 
needs a bounce, as in the Pre Big-Bang model studied in 
the Einstein frame. Unanimous conclusions on the re- 
sulting spectrum of metric fluctuations in the expanding 
stage are still to come. Among the ongoing controversies, 
one concerned the fate of cosmological perturbations in a 
hydrodynamical radiation bounce triggered by a negative 
energy density scalar field 0, generalized afterwards in 
. This bounce - and its generalization - has been sug- 



gested as a simple toy model, which has the advantage of 
providing analytic solutions for the background, although 
having a component which violates the null energy condi- 
tion. Note that without assuming spatial curvature and 
demanding general relativity to hold, such a negative en- 
ergy component is required at the level of an effective 
theory in order for a bounce to take place Q . 



The initial result in this class of two-fluid models was 
that the spectrum of the Newtonian potential after the 
bounce was the same as that of the growing mode be- 
fore the bounce. Such a result was obtained evolving 
numerically 0] and analytically 0, HI a set of regular 
equations. This result was later challenged in [l^, 11 1, 
generalized in Ref. [l^ in which the scalar perturbations 
are evolved through a bounce characterized by a sin- 
gle physical scale, arguing that the growing mode before 
the bounce matched only with the decaying mode after 



the bounce, a possibility which has been already found 
[S H Q El- However, the analysis of [13, [ill, de- 
manding the most general possible situation (the case at 
hand in the present work being a subset), needs to rely on 
a set of singular equations, a fact that could cast doubts 
on its accuracy had they use them directly; these authors 
however obtained the solution in the form of a Born-like 
series containing only convergent integrals; Non-singular 
equations have also been evolved in different contexts, 
e.g. with a double scalar field bounce [lB| or one with a 
non-local dilaton potential stemming from string theory 
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17| . Note that for tensor modes, we already know that 
a matching between growing modes before and after the 
bounce occurs [1, • 

In this paper we consequently reanalyze the be- 
haviour of cosmological perturbations during the radi- 
ation bounce, obtaining results in agreement with these 
later studies 13, Hi, and in contrast with our previous 
findings for scalar perturbations in 0, Q • In section II 
we present the background bouncing models containing 
two perfect fluids with linear and unrelated equations of 
state. We also show how to describe the negative en- 
ergy perfect fluid in terms of a K-essence scalar field. In 
section III we propose a set of regular equations for lin- 
ear perturbations of the above background models, one 
for the Newtonian potential and the other for the ve- 
locity potential of the fluid responsible for the bounce, 
which can be related to the linear perturbations of the 
K-essence scalar field yielding simpler regular equations 
suitable for the numerical analysis presented in section 
IV (another set involving the density contrast instead of 
the velocity potential is given in the appendix). In sec- 
tion V we justify some of the numerical results through 
an analitical study of approximate solutions and their 
matchings. We end in section VI with discussions and 
conclusions. 
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Note that the new coordinate time r makes it possible 
to generalize the solution, obtained in terms of the usual 
conformal time found in [8] to arbitrary values of w+ , w- . 
Note also that this new coordinate time allowed to get 
general solutions for the scale factor in a universe filled 
by dust plus dark energy, the latter having an arbitrary 
constant equation of state [lH . 

One can also describe fluids in terms of velocity po- 
tentials [13]. In the case of a perfect fluid, the velocity 
potential action is very simple, and identical to an action 
for a K-essence scalar field [21| (which, incidentally, can 
also be used to produce a bounce with no curvature and 
only one scalar degree of freedom namely 



S 



-gd X, 



(8) 



BACKGROUND 

We shall consider a class of bouncing universes filled 
by two non-interacting perfect fluids with parameters of 
state 1/;+ , W- , constant in time 8] , relating the energy 
densities e± to the pressures p± through p± = w±e±. 
The Einstein energy constraint for homogeneous and 
isotropic solutions is 



aat , 



P+ 



P- 



(1) 



where SttG — M^^ — S^^^, p+, p_ being constants. Eq. 



([T]) is obtained using the background FLRW metric 



As' 



(2) 



and by assuming energy conservation for each single fluid 
separately in order to make explicit its dependence on the 
scale factor. We restrict ourselves to W- > w+. It is clear 
that the negative energy density fluid, or, in other words, 
p_, is important only close to the bounce, in agreement 
with what is required from a phenomenological model. 
By introducing a new coordinate time t: 



dr ^ -g . 



with 13 = -{2W+-W- + 1), (3) 



we can solve Eq. ([TJ for the scale factor as 



a(T) = ao 



T 

1 + 



(4) 



where the ± sign is chosen according to whether the 
fluid has positive or negative energy density. The energy- 
momentum tensor for 6 reads 



± 



(1 +w±) 



2w± 



'5w1,A(1-'^±)/(2«'±) 



(9) 



and the field equation of motion (stemming from the 
energy-momentum conservation) is given by 



(1 - w±) W^cbW^cbW^Wpcbg^'^g'^f' 



w± 



0. (10) 



In the homogeneous background the energy den- 
sity and pressure read 



e±(o) = T'^oo - ± 



1 



2w± \ a 



I \ {l+w±)/w± 



P±io) 



(0), 



(11) 



(12) 



3 ~^^(o)o 

and the equation of motion JTU]) reduces to 

ip" + {3w± - IW = 0, (13) 

where tp, e±(o) and p±(o) are the homogeneous parts of 
(j)^ e± and p±, respectively. In pTjl and the following 
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equations, a prime represents a derivative with respect 
to the conformal time rj of the metric (fTT)) . 

From the equation for tp one obtains tp' = 
C^/a^^^^~^\ where is a constant related to p± 
through 



In order to make contact with the notation of Ref. [7|, 
which is contained in the general case we treat here (a = 
1/2, w_ = 1), we will choose to represent the negative 
energy perfect fluid by this K-essence scalar field and let 
the positive energy fluid with its original hydrodynamical 
representation, thus leading to the action 



5 = - 



1 



167rG 



R- 



•e+ + -(V^0V 



A'^^(l+™-)/(2'^-) 



(15) 



r 



where R is the curvature scalar that takes into account 
gravity. Eq. ([T5)l reduces to that of Ref. 0] if W- = 1- 
The full Einstein equations then reads 



Gfxu — 87rG'(r+ + T ), 



(16) 



form of scalar metric perturbations on the background 
given by Eq. ^ reads, in the longitudinal gauge, 

ds^ = a^{r^) [(1 + 2$)d772 - (1 - 2$)%da:Mx^] , (17) 



where G^tiy is the Einstein tensor, are the energy- 
momentum tensors of the positive and negative energy 
perfect fluid respectively, with written in the hydro- 
dynamical representation, and T'^^ expressed in terms of 
(/) given by Eq. ([9]) . 



Regular equations for cosmological perturbations. 

A set of regular equations is a necessary tool for a 
numerical analysis of bouncephysics. We shall gener- 
alize the treatment of Ref. [7|] to the generalized class 
of bouncing models found in Sect. II. The most general 



J 



where $ is the gauge invariant Bardeen potential j23l.l24|. 
For the matter fields we have 



= <^('7) + '^0(x,?y) and e± = £±(77) + '^e±(x, t?), (18) 

where e± ee ±T^^. 

From Eq. ([9]) at first order we obtain 



(5_ = 



5^ 

V 



- $ 



(19) 



One can also check from Eq. that (5p_ = ST^/3 



Using Eq. ()10|1 to obtain a linear equation for (50, and the perturbed Einstein equations in order to obtain an 
equation for the Bardeen potential after eliminating 5p+, yields the coupled set of regular equations for the modes 
of wavenumber k, namely 

w+k^ + 2H' + (H^ - K){1 + 3w+) + ^H^Q^F $i 



$'fc'-f3H(l + zi;+)$; 



In'n.F^-^, (20) 
2 ip 



and 



b'l + H(3w_ - l)6(t)'k + k^w^5(t>k = (1 + iw^y^'k , 



(21) 



which we wrote in full generality by including a possibly non- vanishing curvature of the homogeneous spatial section 
K, and we have set F = (w+ — w^){l + w^)/w^, and = e^i'^^a'^ /Ti'^ . 

Another way one can write Eq. (|20p using the background equations of motion, which will be useful! when discussing 
the possible spectra, reads; 



-1-37^(1 + ^;+)$',. 
3 



2w- 



(1 + 3w+) 



(3w- + 1) 
2w- 



K 



4m;_ 



e (2w;_p_)i/(i+»-)F^. 



(22) 



The advantage of Eqs. (|20I21[) is their general use in bounces with two hydrodynamical fluids, constituting a 
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set of coupled regular equations for numerical analysis 
of bounce physics. We should like to use this opportu- 
nity to mention the fact that attempts to write down two 
uncoupled equations for two separate parts of the New- 
tonian potential, as suggested in Ref. 0] are incorrect 
ii, 111,123, as well as its use Q. 

Another set of regular equations for general two-fluid 
models solely in terms of hydrodynamical variables (with- 
out using the K-essence scalar field to describe the 
W- perfect fluid) can be obtained using, together with 
the perturbed Einstein equations, the perturbed energy- 



momentum tensor conservation equations 

6'_, + {l+w^){e'_,-3<P',)^0, (23) 



0L,+7^(l~3^;_)g_fe-fc^$fc^ ^-^"^-fe , (24) 

1 + w_ 

where {6uk-)i = adiO^k/k"^ [(<5ufc-)i is the perturbed 
W- three- velocity mode], and (5_ = 5p^/ p^. The result 
reads 



$',: + 37^(l + u;+)$',. 



w+k'^ + 2n' + {n^ - K){i + 3w+) + -n^n^F 



fc2 



6*" +7i:(l - 3w-)9'_k+ [k^w- + (1 -3w_)7i:'] 0-k = fc^(l + 3w_)$'fe , 
where F = {w+ — W-){1 + W-)/w- and il- = p-^l^ /Ti^, as before. 



(25) 
(26) 



The relation between 9 and Scj) is given by. 



NUMERICAL RESULTS 



(27) 



from which the system (I20l2ip can be recovered from 
(|25l26p straightforwardly. 



The system (|20l2ip written in terms of and Yfc — 
Sipipi which evolve in the variable x ~ t/tq for the family 
of bounces described by ^ read 



dx^ 



+ (1 + a + ap)^ 



«;+fc2(l + a;2)M/3-i) + 



4q!(1 -3a + afi) 
(1 + 2:2)2(2- 3a -H2a/3) 



2\/2a(a/3 + 1)(1 - a + a(3) 
2 - 3a + 2aP ^ 



dec 



for the metric perturbation, and 



da;2 



+ (2 + a/3 -3a) 



.x2 -\- \ Ax 



a 



Ax 



(28) 



(29) 



with k = kTQGQ^^ for the scalar fleld part. 

In what follows below, we have solved, numerically, the set consisting of Eqs. and ([^5]) . setting unnormalized 
vacuum initial conditions (the fact that we do not bother about the normalization here is because we are merely 
interested in the transmitted spectrum), reading 



2,-3q(1+iu-)-) 

*fe,ini = 7= exp 



■ K^W+ ^2a(/3-l)-H 



X 



a(l — S-u;— ) 



Vk 



■ exp 



l-^2a(/?- 1) 
k,/w^ 



l + 2a(/3-l) 
I 



2a(/3-l) + l 



(30) 



Fig. [T] shows the time evolution for the spectrum cal value for the scalar spectral index ng in the expanding 
P$/fc"^~^, for that particular case for which the theoreti- 
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stage is known and given by Eq. (|5T|) . i.e. for w- = 1 
and = 10~^. The plots for three different wavenum- 
bers show how the spectrum for the (growing mode of 
the) Newtonian potential in the expanding stage is the 
one of curvature perturbations in the contracting stage. 
Fig. O shows the same plot (near the bounce only) for a 
different value of W-, namely w- = 1/4, again rescaled 
with the theoretical prediction. 



ANALYTIC STUDY OF THE STIFF MATTER 
CASE 

In this section, we will analytically justify the spectra 
of the class of bounces which are driven by negative en- 
ergy stiff matter. Note that for w_ = 1 in the equation 
for the Newtonian potential $fc (|22|) . the "mass" term 
on the left hand side vanishes for flat spatial sections, 
K = 0, and infinite wavelength, fc = 0. In this case, 
one can obtain a solution around the bounce which can 
be matched with the solutions far from it to obtain the 
spectra for large wavelengths. This section is a general- 
ization of what was done in Ref . [7] , which concentrated 
on the radiation-stiff matter case only. 



The relevant phases in the perturbations evolution 

Let us first consider the asymptotic limit r; — > — oo, 
where x oc r^3(i-™+)/(i+3«)+) j^^^j ^ ^ ^2/(i-h3™+) ^^^le 

positive energy fluid dominates). Taking into account 
the initial conditions pO|) . which in terms of rj read 



from which one derives and 



exp{—ikrj) 
exp(— z/cry) 



(31) 



one can see that the source terms in Eqs. ([22l and ((2T|) 
are negligible provided < 7/3. Defining the variables 
Uk = a3(i+'"+)/2$fc and Vk = a^^'"—^^'^5(j)k, one obtains 
the equations: 



6 



and 



(l + 3w+)V 

,(3™--l)/2-)"" 



j(3«;_-l)/2 



Uk = 0, 



Vk = 0. 



(32) 



(33) 



+ X,,,<'(ViyrA;ry)]. 



(35) 



where v 



5+3u;+ 



and /i 



The coeffi- 



2(l+3to+) ' 2(l+3u)+) 

cients '^^.■^ and X^.^ arc time-independent and only de- 
pend on k. 

Restricting now to the case W- — 1, and taking 
into account the initial conditions (j3ip . one obtains that 
<I>jjj = X^^.^ = 0, = 1/fc, and X^^^ = 1. In the region 
where krj <C 1 (we are considering long wavelengths) but 
still far from the bounce, where the source terms can still 
be neglected, one has 



^<k=Ai 



+ 0(eif), (36) 



and 



9< k 



B2 

^3(l-tu+)/(l+3i«+) 



-HO(fcV), (37) 



where Ai = fc3(i-t.+ )/[2(i+3»+)]^ 

^(-7-9u.+ )/[2(l+3«)+)]^ ^ ^3(l-»+)/[2(l-|-3i«+)] ^-^^ 

B2 = fc3(«'+-l)/[2(l+3'^+)l. 

Now we have to propagate this solution through the 
bounce and match it with the solution in the expanding 
phase. During the bounce, the source terms in Eqs. ([22]) 
and pT|) cannot be neglected but now it is the terms pro- 
portional to k^ that are negligible. For w- = 1, Eq. 
simplifies and, upon returning to the variables y = £p 
and X = t/tq, one obtains 



d^$fe_^ Sax d$fc 
dx^ : 



V2 m 



1 da; (1 4- x^) da; 



d^Yfe 2x dYfe _ 8V2a A<^>k 
da;^ a;^ + 1 da; (1 + x"^) dx 



(38) 



(39) 



where, as W- = 1, a = 1/[3(1 — w+)]. The solutions of 
these equations read 



i + B/i(a;) + C/2(a;), 



(40) 



and 



The solutions of these equations are 

Uk = Vv ^(1) hI'^ (\/w+fcf?) + *(2) -f^r {^/w^kvi) 



Y^°^--- = b + Bh{x) + Cfi{x), 



(41) 



' with A, i?, C and D arbitrary constants. The bounce 

(34) functions fi{x) are found to be 



I I \ I I I I til I I I I I I I I 

-50 50 100 -0.01 -0.005 0.005 0.01 

X X 

FIG. 1: Example of the time dependence of the Newtonian potential for three different wavelengths (fc — 10~^, lO"'^ and 10~^ 
respectively) as function of a; = qjjjo- This example, for which w- — 1 and w+ = 10~^ is typical of most cases for which there 
is a constant mode, as found in [lol . [TH . Note that the amplitude at the bounce can be much larger than that of the constant 
mode that dominates later. The Bardeen potential is here rescaled by the predicted spectrum, which in this case is given by 
Eq. (EIJ 



(42) 



f2ix) 



8a / dx 



(l+i2)l+4a 



^,2 



l)F[-,-4a- 



'2' 



-^2 



— , —4a, — , —X 
2' ' 2' 



(43) 



dx 



(1+^2)1+40^ V 2' ' 2' 



where F denotes the hypergeometric function. For a; 3> 1, the solutions can be written as 



(44) 



A- 



B 



C 



r8a~l 



A 



B 



C 



j^{b+3w+)/(l+3w+) ^6(1— uj+)/(l+3i«+) ■ 



and 



yBounce ^ f> - + £ ^ £) 



B 



C 



j^8/(l+3i«+) j^3(l-u)+)/(l-l-3t«+) ■ 



These solutions coincide with those obtained in Ref. for w+ = 1/3. 



(45) 



(46) 



If we now compare Eqs. (|45l46p with Eqs. ()36|37|) . 
one can obtain that A = /s3(i-«'+)/[2(i+3«.+)] ^ ^ ^ 

^(--7-9u;+)/[2(l+3i«+)]^ f) ^ ^3(l-i«+)/[2(H-3«)+)] j^^^j (j ^ 

k^{w+-i)/[2(i+zw+)]^ One can also see this by noting that 
the third term in Eq. (|45p is the first contribution of 
the source term to $, which of course must have the 
/c-dependence of B2 in (|37|) . while the second term of 
Eq. ([^5]) is the first contribution of the source term to 



which must have the /c-dependence of A2 in (|36|) . 



We now have to match Eqs. (|45|46p with the solutions 
in the expanding phase which are far from the bounce 
and where the source terms are negligible, i.e.: 



> k 



jy(5-f3u)+)/(H-3u.+) 



^OikS^), (47) 
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FIG. 2: Same as FIG. [T]with a different value for 10+, namely 
w+ = 1/4. Again, the Bardeen potential is rescaled with 
the spectrum found by [13, [Ul, which we thus independently 
confirm. 



If we now calculate the power spectrum 



we obtain 



1 = 



12w+ 
1 + 3w+ ' 



(50) 



(51) 



as obtained in the numerical calculations and in Refs. 
ll| . This result was also obtained by considerations on 
the matching conditions in 0, [l3| , which predict the 
spectral index in the expanding stage as the one of cur- 
vature perturbations in the contracting stage. Note inci- 
dentally that it coincides with the spectrum obtained in 
Ref. L28!|, where the bounce is not caused by a negative 
energy stiff matter but by quantum effects: the back- 
ground and the spectrum have the same behavior. 



and 



GRAVITATIONAL WAVES 



Bo 



^3(«;+-l)/(l-H3t«+) 



+ 0(A:V), (48) 



yielding, for the constant part of which determines the 
spectrum, 

Ai OC fc3(l-to+)/[2(l+3u.+ )] _^ ^3(™+-l)/[2(l-H3u>+)] 

^ ^3(«;+-l)/[2(l-H3u.+)]^ (-49^ 

because 10+ < 1 and we assume fc ^ 1 . Note that $ gets 
the spectrum of i5(/). 



The equation for the Fourier transforms of the am- 
plitude of the two polarization degrees of gravitational 
waves in cosmology is 



di2 



dt 



^2 



(52) 



Once we introduce v = a'^^~^^/^h and use the same co- 
ordinate variable r as introduced in Eq. (2), the above 
equation becomes: 



J 



Vk + 



-1) , {I3~3)d (/3-l)(/?-3) fd 



Vk = 0, 



(53) 



where / = dj/dr. The very existence of an asymptotic vacuum demands the condition 2a(l — (3) < 1, or, in 
other words, if w+ > —1/3, which we shall therefore assume. From here on, for convenience, we also define S = 
1 -I- 2a{f3 — 1) > 0. We also restrict to w+ < 1 in order to have the constant mode as the growing mode in the 
expanding stage. 

The equation which we numerically evolve is: 



dx2 



2a(/3-l) _^ a{l3-3) 



Vk = 



(54) 



with, as usual, x = t/tq, k = ka^ ^tq. 

Finally, the initial conditions corresponding to the adiabatic vacuum are taken to be 



^fk 



exp(— ifcry) 



,ini = \/3to^p 



j(l-/3) 



\fk 



■ e'Kp{—ikri), 



(55) 
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where = a hk, and we get rid of the prefactor since we are mostly interested in the spectral index anyway (just like 
for the scalar case, the normalisation here is essentially irrelevant). Recall also that kr] = kx^ /S. 



Analytic Approximations 



We are first interested in determining the matching point between the short and long wavelength approximations. 
The potential in terms of the conformal time is 



n" 

a 



- + (1-/3) - 
a \a 



Expliciting this in the r variable, this is 



where 



V = 2aaf -''V(r), 



fir) . - {[1 + 2a (/? - 2)] r„^} rt'^-^^ (r^ + ,2) 
so the matching point at which k^ ^ \a" /a\, i.e 



" [ V2« [1 + 2a(/3 - 2)] j 

where the last inequality stems from the requirement that there is an asymptotic vacuum, i.e. S > 0. 
The zeros of the first derivative of V are determined by the equation 

T{[2a(4a - 3 - 6a/3 + 2^ + 2a(3^) + 1]t^ + [2a{3 - 2/3) - 3]t^} = 0. 

I 



(56) 



(57) 



(58) 



(59) 



(60) 



We will here treat the simplest case where the potential 
V has only one extremal point, at r = 0, hence imposing 
that the coefficients of and Tq have the same sign. 
Asymptotically far from the bounce, Eq. ([54|) becomes: 



dx^ 



■p^4c.(/3-l) _ 7(1 + 7) 



Vk = 0, 



(61) 



where 7 = a(/3 — 3); the above equation admits a solution 
in terms of the Hankcl function, in accordance with the 
vacuum initial conditions (15511: 



(2) 



(62) 



supposed to be valid up to a;„ of Eq. (|59p. with 

A^ = ^^^,exp-('^+i), 



and 



7 + 



Note incidentally at this point that the matching time 
(I59p gives an argument for the Hankel function which 



does not depend on k. We shall henceforth call 
Hi^\~kxl/5)=Hi^\l/6)^K'. 

On the other hand, for long wavelengths close to the 
bounce, Eq. ([51)1 simplifies to 



d^Vk , 7[i-(i + 7); 



da;^ 



{1 + x^y 



-Vk = 0. 



(63) 



In this limit, setting Vk = ^/l-i-x^u and z = ix, one gets 
the Legendre equation 



d^u 

dl2 



-2z 



dz 



7(7 + 1)- 



(1+7)' 



u = 0, (64) 



which in this case has, as the two independent solutions, a 
power law and a hypergeometric function. Summarizing, 
we obtain, in this second regime, the general solution 



Vk,2 



1 + x' 



-7/2 



B + CX2F1 { 



B + C 



V^r(-i-7) 
2r(-7) 



1 + 27 



(65) 
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so we end up with 



FIG. 3: Example of the time dependence of the tensor per- 
turbation for three different wavelengths (k = 10~^, 10~^ and 
10^^ respectively) as function of ki]. This example, for which 
w- — 1/3 and w+ — 10~^ is typical, again, and in complete 
agreement with our analytical prediction. All other cases lead 
to identical figures, except for the actual numbers. The cases 
for which the potential does not satisfy the condition that 
the potential has only one extremum at r = cannot be 
compared with theoretical expectation, in the lack of it, and 
so are not shown. 



It is now a simple matter to match the solutions 
and (|65p as well as their derivatives to get 



Ah^'k-^/'^^^^ = {B + CT)k^/^ + 



C 



1 + 27 



^(-i-7)/5^ (66) 



and 



= (-CT - 5)7^(^+1)/'^ + — ^^-'^/^ (67) 



1 + 27 



where we have set 



V^r(-i-7) 
- 2r(-7) 

for notational simplicity. The solution of this system 
provides B and C a functions of the reduced wavenum- 
ber /c, and we shall retain in what follows the leading 
order terms, which is, as we are considering wj^ < 1, 

^(1+27)7(25) ^ p{w+~l)/[2i3w+ + l)]^ yielding for h W X^^i, 
h « ^3(™+-l)/[2(3^+ + l)] (^Qjjg^^ ^ x^-^ + l). (68) 

The actual gravitational wave spectrum is 



being 



12i«4 



2a 



l + 3w+ l + 2a(/3-l) 



(70) 



(71) 



It is worth pointing out at this stage that Eq. (|7ip gives 
the same result as in the scalar case [Eq. ([51]) ] for the 
specific case that we could study analytically. The rea- 
son for such similar results stems from the fact that the 
dominant terms which match through the bounces under 
investigation are the growing modes of curvature per- 
turbation and gravitational waves, both satisfying the 
same differential equation in the single fluid regime. The 
above result (I7ip was already obtained in previous inves- 
tigations 8, 28], although for two different subsets of the 
family of bounces studied here. 



CONCLUSIONS 

In all the early universe models which aim at solv- 
ing the horizon problem with a contraction instead of 
a superluminal expansion, a deep understanding of the 
physics at the bounce is crucial (and presently lacking 
in its full generality). What we have shown here is a 
step towards the understanding of cosmological pertur- 
bations through a bounce triggered by a second perfect 
fluid (with negative energy density), in the framework of 
flat spatial section and general relativity. 

We have analysed in greater details, both numeri- 
cally and analytically, this class of two perfect fluid 
bounces with flat spatial sections, using a completly reg- 
ular system of equations, concluding indeed that the con- 
stant mode of the scalar gravitational potential after the 
bounce does not acquire a piece of the growing mode be- 
fore the bounce. Therefore, our conclusions agree with 
ill ] , and are in contrast _with our previous results 



\h\\ 



(69) 



[1 ^, 

for scalar perturbations in 0, [1] ■ Another important re- 
sult is the unsensitivity of the scalar spectral index from 
the peculiarities of the bouncing component in the class 
of models studied in this paper. One interesting result 
is that when the negative energy fluid has stiff matter 
equation of state, the background model and the pertur- 
bations have the same behaviour as the quantum bounc- 
ing cosmological models analyzed in Ref. J8]. Our re- 
sults are interesting for the predictions of cosmological 
alternative models. Whereas by a very slow contraction 
- as in Ekpyrotic/cyclic model - it seems really difficult to 
generate a nearly scale-invariant spectrum of curvature 
perturbations without the need of isocurvature perturba- 
tions or extra-dimensions, a homogeneous dust contrac- 
tion 13, [28| seems in agreement with observations and 
even free from details due to the bouncing component 
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which were left open from previous investigations which 
focuses on w- = 1 [3 28 1. Note however, as far as 
complete model building is concerned, the assumption of 
homogeneity may not hold close to the bounce and should 
thus be verified. This point is however out of the scope of 
the present article whose aim was to concentrate on the 
propagation of scalar and tensor perturbations through 
a regular, although phenomenological, bounce. 

We have also performed the analytical and numerical 
calculations for gravitational waves. In this case, the 
constant mode of the long wavelength tensor perturba- 
tions after the bounce do acquire a piece of the growing 
mode before the bounce. Also in this case, the slope of 
the final spectrum does not depend on the negative en- 
ergy perfect fluid equation of state. This paves the way 
to a generic behavior for tensor perturbations, as such 
a phenomenological model thus does not suffer from the 
drawback (still present for the scalar modes) of relying 
heavily upon the details of the bounce physics. Both 
these results agree with the previous investigation Q for 
a restricted class of bounces. This can be understood by 
noticing that the crucial time in the evolution of the per- 
turbations is when the perturbation wavelength becomes 
comparable with the curvature scale of the background, 
when, for large wavelengths, the universe is still far from 
the bounce and hence the effects of the negative energy 
fluid are negligible. This result was already anticipated 
in Ref. 
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APPENDIX: REGULAR EQUATIONS FOR 
GENERAL TWO-FLUID MODELS IN TERMS OF 
HYDRODYNAMICAL VARIABLES. 



Another possible set of regular equations uses the den- 
sity contrast 5- = Sp- / p_ of the fluid driving the bounce 
instead of its velocity potential. The equations are: 



3H{l + w+)^'k + {w+k^ + 2H' + + 3w+H^) <i>k = -7i^(5_0_(w+ - w_) 



(72) 



and 



5'1 +{l-3w-)n6'_ 



il + w^){ [fc2(l + 3w+) + 3[2n' + (1 + 3w+)n^] $fe + 3[2 -I- 3( 



■ w. 



(73) 



We note that with this new set the order of the system of linear differential equations is increased with respect to the 
systems (|20l2ip or (|22l26p : Eq. ([73]) is indeed equivalent to a third order differential equation for Scfik [see Eq. ([T9)) ]. 
As a result, solving this last set of equations may lead to spurious solutions and it is therefore better to stick with 
Eqs. ((20|2T|) . 
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